Magnetic field induced phase transitions in spin ladders with ferromagnetic legs 
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We study the phase diagram of spin ladders with ferromagnetic legs under the influence of a 
symmetry breaking magnetic field in the weak coupling effective field theory by bosonization. For 
antiferromagnetic interleg coupling we identify two phase transitions introduced by the external 
magnetic field. In order to establish the universality of the phases we supplement the bosonization 
approach by results from a strong coupling (rung) expansion and from spin wave analysis. 
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I. INTRODUCTION 

Recently there has been considerable interest in the 
study of magnetic field-induced effects in low-dimensional 
quantum spin systems, in particular devoted to the 
critical properties of spin S = 1/2 isotropic antiferro- 
magnetic two-leg ladders in an external magnetic field. 
In parallel, there was remarkable progress in recent 
years in the fabrication of such ladder compounds^. 
Since antiferromagnetic two-leg ladder systems with 
S = 1/2 have a gap in the spin excitation spectrum, 
they reveal an extremely rich behavior, dominated by 
quantum effects, in the presence of a magnetic field. 
These quantum phase transitions were intensively in- 
vestigated both theoretically, using different analytical 
and numerical t echnique a^ • ^ S ' i " ! "*^ i -'^ ^ • ""^ ? i "'^ i "'^ ^ i "'^ ? , and also 
exp er iment all jiSiMiSiSil. 

Ladder models with ferromagnetic legs have been much 
less studied, although they exhibit many interesting 
aspectsiiiiSiiS. It is true that up to now no materials 
are available which realize these models. However, from 
the theoretical point of view these systems are extremely 
interesting, since they open up a new large class of sys- 
tems for the study of complicated quantum behavior, un- 
expected in more conventional spin systems. The vari- 
ety of possibilities is seen already from Fig. ^ here the 
ground state phase diagramiS of a two-leg ferromagnetic 
ladder is presented in the variables intraleg exchange 
anisotropy (A) and (isotropic) interleg coupling (Jx). 
The ground state phase diagram contains, besides the 
fully gapped rung-singlet and Haldane phases (commonly 
known from the case of antiferromagnetic laddersi^-'^°), 
the spin-liquid phase with easy-plane anisotropy (XKl), 
the ferromagnetic and the stripe-ferromagnetic phases 
which are realized only in the case of ferromagnetic legs 
(0 < A < 1). 

In this paper we study the effect of an external mag- 
netic field on the phase diagram of this system. In par- 
ticular, we focus our attention on the study of new field 
induced effects in the case of the easy-plane XYl phase 
and in the rung-singlet phase in the vicinity of the single 
chain ferromagnetic instability point A = 1. 

The effect of the uniform magnetic field, applied paral- 
lel to the anisotropy (z) axis, on the ground state prop- 



erties of the two-leg ladder systems is known from the 
investigation by Schulali for the case of the correspond- 
ing spin 5 = 1 Heisenberg chain model: In the case of 
the gapless XYl phase an external magnetic field leads 
to the appearance of a magnetization in z direction for 
arbitrary small magnetic fields. In the case of the gapped 
rung-singlet phase the magnetization appears only at a 
finite critical value of the magnetic field which is equal 
to the spin gap^"^^. This behavior is generic for the spin 
gapped U{1) symmetric systems in a magnetic field: The 
magnetic field leaves the in-plane rotational invariance 
unchanged?! and the transition belongs to the univer- 
sality class of the commensurate-incommensurate (C-IC) 
transitionsSSiS^. 

On the other side, for the U{1) symmetric phase such 
as the XYl phase, the effect of a uniform transverse field 
is highly nontrivial. In the case of classical anisotropic 
spin chains this effect has been studied more than two 
decades agc^^. However, in the case of the antiferro- 
magnetic XXZ quantum chain this problem is still the 
subject of intensive recent studie32^*2Si2L2&. 

In this paper, we study the effect of a uniform trans- 
verse magnetic field on the ground state phase diagram 
of a two-leg ladder with anisotropic, ferromagnetically 
interacting legs coupled by antiferromagnetic interleg ex- 
change. 

The outline of the paper is as follows: In section II 
we review the model and its bosonized version in the 
continuum limit. We discuss the phases and phase tran- 
sitions emerging from the XYl and rung-singlet phases 
in a transverse magnetic field in sections III (in the weak 
coupling approach) and IV (in the limit of strong rung 
exchange). We shortly summarize our results in section 
V. In appendix A we present the spin-wave approach to 
locate the ferromagnetic transition line starting from the 
saturated phase. 



II. MODEL 

Here we present a brief introduction to the model and 
in particular to its bosonized version in the contimmm 
limit. The Hamiltonian we consider is given by: 



(1) 

leg 



(2) 
leg 



(1) 



2 



FERROMAGNET 



A A "Stripe" 
FM 




RUNG 
SINGLETS 



Anisotropic Haldane 
Phase 



RUNG 
SINGLETS 



'1 



FIG. 1: Schematic picture of the ground state phase dia- 
gram of the two-leg ladder in the variables intraleg exchange 
anisotropy (A) and isotropic interleg coupling (J_i). The 
limiting case of isotropic ferromagnetic legs corresponds to 
A = 1. 



where the Hamiltonian for leg a is 
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and the interleg coupling is given by 
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Here S^''^' are spin S = 1/2 operators on the j-th rung, 
and the index a = 1,2 denotes the ladder legs. The 
intraleg coupling constant is ferromagnetic, J > 0, and 
therefore the limiting case of isotropic ferromagnetic legs 
corresponds to A = 1. We will restrict ourselves to the 
case < A < 1. 

We use the following bosonization expressions for spin 
operators, adapted to the case of ferromagnetic exchange 
(for details seei*^): 
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Here, (j){x) and 9{x) are dual bosonic fields, dt(j) = ud^d 
and K is the Luttinger liquid parameter 
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Now we introduce the symmetric and antisymmetric 
combinations of the bosonic fields (/)± = ^1/2 {ipi ± (^2), 
d± = Y^l/2 (9i ± 62) and after rescaling these fields we 
obtain finally as effective bosonic Hamiltonian 



n = H+ + H- +W 



(8) 



where 



n+ = 



^[(dxO, 



cos ^/8^^K^(j)-^- (x) , 
— [id,9^)^ + id,^_ix)f] 
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J±cos . —e_{x), 
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Here 



Jj_ 2K -1 \ 



in/2K)) ' ^^^^ 



27r J sin 



J7j_ ^ J±, h ^ /i*^^* and u± are the velocities of the 
symmetric and antisymmetric modes. 

In deriving (jSJ, several terms which are strongly ir- 
relevant in our case of a ladder with ferromagnetic legs 
and applied transverse magnetic field, were omitted. For 
details of the full Hamiltonian we refer the reader toiS. 

We note that at /i'^^* = the effective theory 
of the original ladder model is given by two decou- 
pled quantum sine-Gordon models which describe, re- 
spectively, the symmetric and antisymmetric degrees 
of freedom. The infrared properties of the antisym- 
metric field are governed by the strongly relevant op- 
erator J7+ cos \/2'K jK-Q- with the scaling dimension 
(2KJ)^^ < 1/2. Therefore, the antisymmetric sector 
is gapped at arbitrary J± ^ 0. Fluctuations of the field 
0-{x) are completely suppressed and the field 0_ gets 
ordered with expectation values 







at Jj_ > 
at J'^ < 
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The infrared properties of the symmetric field are gov- 
erned by the marginal operator M+ cos y/SiTK^cj)^ . As 
shown the symmetric mode remains gapless for fer- 
romagnetic interleg exchange and arbitrary < A < 1, 
while in the case of antiferromagnetic interleg exchange 
it is gapless in a finite regime in A, Jj_ parameter space 
given approximately by 71 Jj_/J < A < 1 — 72J+/J. 
Here 71 and 72 are positive constants (more rigorously 
smooth functions of the anisotropy) of the order of unity. 
Following Schulaii who has discussed a similar phase in 
the context of the spin 5* = 1 chain we denoted this 
phase as spin liquid XYl phase for the spin ladder. At 
J± > 0, outside of this regime, the symmetric mode is 
also gapped and the field 0+ gets ordered and pinned in 



3 



one of its possible minima. We denoted the fully gapped 
phase, realized in the case of antiferromagnetic interleg 
exchange as rung-singlet phase. In this phase spins on 
the same rung tend to form a singlet and the ground 
state corresponds to the state with a singlet pair on each 
rung. The ground state phase diagram of the system at 
f^ext = is given schematically in Fig. ^ 

Below we study the ground state phase diagram of the 
two-leg ladder with anisotropic ferromagnetic legs in the 
presence of a magnetic field. 

III. THE SPIN LADDER IN THE PRESENCE 
OF AN IN-PLANE MAGNETIC FIELD 

In this Section we consider the effect of a transverse 
magnetic field on the ground state phase diagram of the 
spin-liquid XYl phase. Since the magnetic field breaks 
the in-plane rotational symmetry of the XYl phase, it 
is clear that it will eliminate the gapless XYl phase and 
that the system will acquire a gap in the excitation spec- 
trum for arbitrary small h'^^* . Both fields (symmetric and 
antisymmetric) will be pinned in their respective minima. 
In this case the very first direct approach to study the 
ground state phase diagram of the model is to use the 
quasiclassical Ginzburg-Landau type analysis. This ap- 
proach does not cover the limiting case A = 1 where 
bosonization docs not apply. The A = 1 case will be 
discussed in section IV. 



A. Quasiclassical Ginzburg-Landau analysis 

In the presence of an in-plane magnetic field the effec- 
tive Hamiltonian © reduces to two free Gaussian fields 
coupled by the following effective potential 

Ueff (e_,e+) = jj_cos2e_ - /icose_cose+ , (i4) 

where 9± = 7r/2K±d±. To find the vacuum expecta- 
tion values of the pinned field we search for the minima of 
the effective potential Ueff with respect to 0_ and 0+. 
The straightforward analysis gives the following sets of 
vacua: 

u h> hc = 4 Jj_ 

I. = TT , e_ = TT , (mod 27r) 

II. e+ = 0,e_=0, (mod 27r) (15) 

m h < hc = 4:J^ 

I. e+ = TT e_ = TT ± t9o , (mod 27r) 
II. e+ = 0, e_ = ±i?o , (mod 27r) (16) 

where 

•do = arccos(/i/4Jj_). 




FIG. 2: Evolution of the set of minima of the antisymmetric 
field with increase of the applied transverse magnetic field. 

At /i = only the antisymmetric field is gapped and 
its set of available vacua is given by 0_ = 7r/2, (mod tt). 
At ft. 7^ the symmetric field gets pinned and the set 
of possible vacua of the symmetric field does not change 
with field. On the other hand, at arbitrary < h < 4j7j_, 
each minimum in the antisymmetric sector splits into two 
degenerate minima (see Fig. At /i — > he, i^o '^/^ 
the split minima join each other and form a new set of 
possible vacua for the antisymmetric field which is fixed 
for arbitrary h > he- At the critical point the effective 
potential transforms into the (O-)"* potential, which is 
common in describing the Ising universality class. Since 
the location of the minima in the symmetric sector does 
not change, we conclude that the transition involves only 
the antisymmetric sector. 

An analogous quasiclassical analysis was carried out 
by Fabrizio et al.^^ in connection with the ionic Hub- 
bard model. There two transitions were identified for 
finite values of alternating on-site energies. We note as 
an important difference that assuming a vanishing order 
parameter along the direction of an infinitcsimally small 
applied magnetic field would be self-contradictory in our 
case, because here no appropriate marginal operators are 
present. 

Note that in the case of ferromagnetic interleg ex- 
change (J_L < 0), the set of minima of the effective poten- 
tial for arbitrary h^^* is given by the set (|15l) . There- 
fore, no transitions with increasing field are expected in 
this case. 



B. Phase diagram 

The quasi-classical analysis performed above allows to 
sketch qualitatively the ground state phase diagram of 
the model under consideration. 

As soon as we switch on an infinitesimally small mag- 
netic field in X direction, breaking the in-plane rotation 
symmetry, the XYl phase is expelled, the system be- 
comes gapped and acquires a finite order parameter. Us- 
ing the bosonized expressions for the smooth parts of 
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FIG. 3: Crossover from stripe-ferromagnetic (a) to ferromag- 
netic (b) states. 

the spin operators Eqs. expressed in terms of the 

symmetric and antisymmetric fields and the vacuum ex- 
pectation values of the corresponding fields ((T3|l - l(TC|) one 
easily obtains that at < /i < ft.c the spin ladder acquires 

• uniform magnetization in the direction of applied 
field 



~ (cose+) (cose_) = V4J± 



(17) 



• opposite magnetization of legs in the in-plane di- 
rection perpendicular to the field: 



~ (cose+) (sine_) = {h/AJ^f . (18) 

This phase we denote as the " stripe- ferromagnetic" 
phase. 

When the magnetic field exceeds the critical value h > 
he, the new set of vacua ((TH|l is reached and the system 
passes into the ferromagnetic phase, where 



(cose_) = land M„ 



(sine) = 0. 



(19) 



Thus, the quasiclassical analysis of the ground state 
phase diagram of the XYl phase in the presence of a 
transverse magnetic field shows a phase transition from 
the stripe-ferromagnetic to the ferromagnetic phase (we 
denote by 'ferromagnetic phase' the phase with the mag- 
netization parallel to the external field as only nonvanish- 
ing order parameter). This analysis also indicates, that 
the transition belongs to the Ising universality class. 

From the above it is also clear that in the case of ferro- 
magnetic interleg exchange {J'± < 0) the ferromagnetic 
phase is realized for arbitrary h ^ 0: This is so, since in 
this case the vacua of two terms in (|14f) do not exclude 
each other. 



C. Double sine-Gordon analysis 

According to the quasiclassical analysis the magnetic 
field induced transition in the XYl phase from the partly 
polarized stripe-ferromagnetic phase into the ferromag- 
netic phase takes place along the critical line he = 4j7j_ 
and belongs to the Ising universality class. In this subsec- 
tion we investigate the stability against quantum fluctu- 
ations of the results obtained so far. For this purpose we 
decouple the initial theory of interacting quantum multi- 
frequency sine-Gordon fields (|HJl to get 



n+ = ^p.0+)2-i-(a.0+)2] 

— cos/3+0+(x) , 

n- = ^p.0_)2 + (9x0-(^))'] 



(20) 
(21) 



- /i_cos/3_6l_.(a;) + J'iCOs2/3_6l_(a;) , (22) 



where 



and 



h± — h{cos \/ t:/2Kj, 



(23) 



(24) 



Thus the symmetric sector is described by the ordinary 
sine-Gordon theory with < 7r/2 for < A < 1, and 
therefore with a strongly relevant massive term. On the 
other hand, in the antisymmetric sector we arrive at the 
double frequency sine-Gordon model with < 7r/2 and 
therefore strongly relevant basic and double-field opera- 
tors. 

In the double frequency sine-Gordon model the quan- 
tum phase transition in the ground state takes place when 
the vacuum configurations of the two cosines compete 
with each other^*^, corresponding to the crossover from 
double well to single well potential in quasiclassical anal- 
ysis. One easily verifies that this criterion can be ap- 
plied analogously to the case of antiferromagnetic rung 
exchange {J± > 0). The dimensional arguments based 
on equating physical masses produced by the two cosine 
terms separately is usually used to define the critical line: 



(h) 
J 



l/{2-dim[h]) 
l/{2-dim[J±]) 



(25) 



Equating these two masses we obtain the following ex- 
pression for the critical line external magnetic field vs. 
interchain coupling: 



her = vJ' 



(26) 



where 



dim[h] 



2 — dim[J±] 



{8K^ - K)/2K+{AK_ - 1) (27) 
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FIG. 4: Appearance of an additional minimum (dashed line) 
which could change the second order phase transition to first 
order, depending on the sign of the higher order harmonic 
generated. 



and r] is some numerical constant of the order of unity. 
The self consistency of the mean field separation insures 
that l(?7|l follows from equating masses produced by mag- 
netic field and interchain couplings separately, before the 
mean field separation. From (|27|l we see her — J± in 
the limit if — > oo, i.e. for the single chain ferromagnetic 
instability point. This will be seen to be consistent with 
the large rung coupling as well as with spin wave analysis 
(see respectively eqs. and (jASIl '). 

Finally we want to mention that the operator prod- 
uct expansion of the two lowest frequencies of the sine- 
Gordon theory for the range of anisotropy parameter 
< A < 1 does not close, i.e. higher relevant harmonics 
are generated in the RG procedure. In some situations 
it is known that higher harmonics can destabilize second 
order phase transitions and make it weakly first ordes^l. 
This happens when generated harmonics introduce new 
minima (e.g. the dashed minimum in Fig. Q) which 
at the phase transition point coexists with the minima 
before the transition. 

For example in the case of only the third harmonic 
retained, the position of the new minimum depends on 
the sign of the generated harmonics. We have followed 
the sign of the generated harmonic along the one loop RG 
flows and checked that in our case there is no indication 
of an upsetting of the second order phase transition. The 
new minima get pinned in between the minima of the first 
two harmonics; they thus may shift the critical line H26|l 
to larger values of the magnetic field but otherwise do 
not introduce any qualitatively new effects. 



IV. LARGE RUNG COUPLING RESULTS 

In this section we consider the effect of an applied mag- 
netic field on the ground state phase diagram of the two- 
leg ladder system with ferromagnetic legs ^ in the lim- 
iting case of strong rung coupling J± ^ J. In this limit 
it is convenient to discuss the model by representing the 
site-spin algebra in terms of on-bond-spin operators'^^. In 
particular, in the case of isotropic interleg exchange the 
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FIG. 5: Phase diagram of a spin ladder with the SU(2) sym- 
metric ferromagnetic legs in a uniform magnetic field. 



ladder can be mapped onto a single S — 1/2 chain^iSiS. 
For completeness we briefly discuss the mapping here: A 
given rung may be in the singlet or in the triplet state 
with energies given by 
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At h'^^^ < J±, one component of the triplet becomes 
closer to the singlet ground state such that for a suffi- 
ciently strong magnetic field we have a situation where 
the singlet and the Sz — +1 component of the triplet 
form a new effective spin 1/2 system. One can easily 
project the original ladder Hamiltonian |^ on the new 
singlet-triplet subspace 



I ^> ^ |t+> = I TT) 

14) ^ |.^> = ^[in>-UT>] 



(28) 



This leads to the definition of the effective spin 1/2 op- 
erators 



1 



'^ri',a=l,2 — ( 1) _y2^" 



(29) 
(30) 



When expressed in terms of the effective spin operators 
(P^ - (|5n|) . the original Hamiltonian (Q) becomes 



H, 



-■'EI? 



I 

—heff -I- Constant, 

i 

where the effective magnetic field 

^e// = h"""* - J_L - - J. 



^ yij , \ X X 1 



(31) 



(32) 



Note that in writing (|31|l for convenience we have ex- 
changed X and Z axis in effective spin space. 

The Hamiltonian (|31|l describes a spin 1/2 fully 
anisotropic XYZ chain in an effective magnetic field. 
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There exists a special case, A = 1, which allows for rig- 
orous analysis. In this case the effective problem reduces 
to the theory of the XX Z chain with a fixed ferromag- 
netic XY anisotropy of 1/2 in an effective magnetic field 
heff. The gapped phase at h^ff < h^}/ for the ladder 
corresponds to the fully saturated magnetization phase 
for the effective spin chain pointing in the direction op- 
posite to the applied field, whereas the massless phase 
for the ladder corresponds to the finite magnetization 
phase of the effective spin- 1/2 chain^. The critical field 
/ig^^ where the ladder is fully magnetized corresponds 
to the fully magnetized phase of the effective spin chain 
pointing along the direction of applied field. From the ex- 
act ground state phase diagram of the anisotropic XX Z 
chain in a magnetic field'^^ using (|32|l we get that the 
isotropic ferromagnetic ladder in a magnetic field shows 
two second order phase transitions: at 

/i""*ci = J± - J (33) 

a transition occurs from the rung dimer to a Luttinger 
liquid phase and at 

h'^'c2 = J± (34) 

a transition from a Luttinger liquid phase into the fully 
polarized phase. 

The transition from the rung-singlet phase into the 
Luttinger liquid phase in the case of the isotropic antifer- 
romagnetic ladder was studied in detail in several recent 
publication3i2iAi*iSi^. It was shown that in the case of 
the gapped rung-singlet phase the magnetization appears 
only at a finite critical value of the magnetic field equal to 
the spin gap. Since this behavior is generic for isotropic 
systems with spin gapSi, and the gap in the ladder sys- 
tem is governed by J± we conclude that the rung-singlet 
to Luttinger liquid phase transition line smoothly reaches 
zero at (see Fig|Sll. 

Note that this large rung coupling analysis reveals that 
the phase transitions in the antiferromagnetically cou- 
pled ladder with ferromagnetic legs in uniform magnetic 
field are connected with those in a ladder with antiferro- 
magnetic legs but in staggered magnetic fieldi£. In both 
cases the magnetic field tries to promote triplets on rungs, 
while the antiferromagnetically coupled ladder supports 
on-rung singlets. 

Away from the isotropic point A = 1 the effective 
Hamiltonian 131(1 describes the fully anisotropic ferro- 
magnetic XYZ chain in a magnetic field that is directed 
perpendicular to the easy axes. For the particular value 
of magnetic field h^ff = 0, the effective XYZ chain is 
long range ordered in direction, corresponding to the 
original ladder system being ordered in the direction per- 
pendicular to the applied magnetic field with opposite 
magnetization on the legs (stripe- ferromagnetic phase). 
For larger values of the effective field it is clear that this 
striped ferromagnetic order will be replaced either by the 
rung singlet phase or the phase with only one order pa- 
rameter - magnetization along the applied field. Thus 




FIG. 6: Phase diagram of a ladder with anisotropic ferromag- 
netic legs in transverse magnetic field. 



we obtain the result that the phases appearing in large 
rung coupling and in weak rung coupling are identical. 
On general grounds one expects an Ising transition to 
take place at critical strengths of the magnetic field also 
for large couplings. Thus for A ^ 1 the Luttinger liq- 
uid phase will be replaced by the striped ferromagnetic 
phase (as in the weak coupling limit), and the transi- 
tions become of Ising type due to the reduced symmetry 
(Fig.©. 

In order to determine the nature of the transition from 
rung singlets to the striped ferromagnetic state in the 
weak coupling limit we can use the fact that in the rung 
singlet phase the operator 



J± cos y'STTK+(j)+ (x) 

is relevant and the ground state consists of nonmagnetic 
singlets situated along the rungs of the ladder. On the 
other hand, while the in-plane magnetic field couples to 
the dual fields (disorder operators) we expect an Ising 
phase transition to take place with the appearance of the 
magnetization perpendicular to the applied field (stripe- 
FM) as dominant order parameter. 



V. CONCLUSIONS 

We have investigated the phase diagram of the 5* = 1/2 
ladder with ferromagnetic legs under the influence of a 
uniform magnetic field breaking the in-plane rotational 
symmetry. In case of antiferromagnetic coupling between 
legs we identified two phase transitions in the plane of 
magnetic field vs interchain coupling. We have extended 
our analysis to the strong rung coupling limit and have 
identified a Luttinger liquid phase which replaces the 
stripe- ferromagnetic phase in the case of SU{2) symmet- 
ric legs. The phase transition line in the case of SU (2) 
ferromagnetic legs was confirmed also by the spin wave 
calculation starting from the saturated limit. 
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APPENDIX A: FERROMAGNETIC 
INSTABILITY 

In this appendix we use the spin wave approach to de- 
termine the critical line corresponding to the ferromag- 
netic instability for the case of SU (2) symmetric legs. We 
refer the reader to the appendix of fl9|) where we have 
considered in detail the analogous analysis of a ladder 



without magnetic field. In the case of SU (2) symmetric 
legs it is straightforward to add the magnetic field term, 
since it couples to the diagonal operator. In this case the 
two sets of spin wave excitation frequencies^^ read: 

u;-{q) = Jcosg-^ (Al) 

cj+{q) = h--t^-I±- J co^ q+-I±. (A2) 

For > we have 

and from the instability condition lu^ (q = 0) = we 
obtain 

/i-^^* = J_L (A3) 
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